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We have studied the propagation of a crack front along the heterogeneous weak plane of a trans- 
parent PMMA block using two different loading conditions: imposed constant velocity and creep 
relaxation. We have focused on the intermittent local dynamics of the fracture front, for a wide range 
of average crack front propagation velocities spanning over four decades. We computed the local 
velocity fluctuations along the fracture front. Two regimes are emphasized: a de-pinning regime of 
high velocity clusters defined as avalanches and a pinning regime of very low velocity creeping lines. 
The scaling properties of the avalanches and pinning lines (size and spatial extent) are found to be 
independent of the loading conditions and of the average crack front velocity. The distribution of 
local fluctuations of the crack front velocity are related to the observed avalanche size distribution. 
Space-time correlations of the local velocities show a simple diffusion growth behaviour. 

PACS numbers: 62.20.mt, 46.50.+a, 68.35.Ct 



I. INTRODUCTION 

Failure of heterogeneous materials has a vast impor- 
tance in geophysical systems, industrial applications and 
of course fundamental physics. This subject is far from 
understood, and has been studied extensively over the 
years [lHH- Of key importance for brittle materials is 
the competition between pinning forces due to local ma- 
terial heterogeneities and elastic forces due to outer ap- 
plied stress, resulting in a complex roughening of fracture 
surfaces. In general this competition triggers a rich his- 
tory dependence of the fracture process. Up until quite 
recently, a broad range of experimental and simulation 
studies have been concerned with the morphology of ei- 
ther fracture surfaces in the case of three-dimensional 
solids [H , or interfacial crack fronts for planar fracture 0- 
7]. In both geometries it has been well established that 
the fracture roughness exhibits self affine scaling prop- 
erties d-QjJ. To this end, theoretical approaches have 
been suggested: the fluctuating line model [III, [IH, where 
the interface is seen as an elastic string propagating in a 
rough morphology, being pinned with different strengths 
at different positio ns, and also the stress weighted per- 
colation approach [lj] with a damage zone ahead of the 
crack. 

In this study, we will pay our attention to the dy- 
namics of fracture propagation. Owing to the mate- 
rial heterogeneities, the motion is complex and charac- 
terised by abrupt jumps separated by periods of rest. 
Both the jumping and the resting behaviour span a large 
range of time scales. This dynamics is often referred 
to as Crackling Noise 1161 . Apart from direct obser- 
vation of fracture |l7J-|2(|, such intermittent dynamics 
embody also large scale activity in earthquakes [2l|- | 23 | . 
acoustic emission during material failure (fiberglass |24| . 



rocks |25| . paper [2q etc.), magnetic domain wall motion 
(Barkhaussen noise) |27|. wetting contact line motion on 
a disordered substrate |28l. [29| , and imbibition fronts in 
porous media [30j . 

Studies on fracture propagation often characterize 
the complex dynamics through related effective aver- 
age quantity, due to the difficulties of direct observation 
and/or insufficient resolution of the spatio-temporal be- 
haviour at local scale. In contrast we use here a trans- 
parent PMMA model for in-plane mode-I fracture well 
suitable for capturing optically detailed intermittent be- 
haviour with high precision in both time and space [f| . 

The present work is a completion and substantial ex- 
tension of the experimental study presented by Mal0y et 
al. in [l8j , where the concept of the waiting time matrix 
was introduced; a consistent way of obtaining the local 
velocity field of the propagation of a pinned interface. 
Statistical analysis, based on the waiting time matrix, 
of avalanche behaviour in fracture front propagation has 
since been followed up by simulations. Bonamy et al. (32j 
quantitatively reproduced the intermittent crackling dy- 
namics observed in experiments, using a crack line model 
based on linear elastic fracture mechanics extended to 
disordered materials. Using a similar string model, but 
with pure quasistatic driving and zero average propaga- 
tion velocity, Laurson et al. [33j have recently proposed 
a scaling relation connecting the global activity with the 
observed local avalanches, connecting the dynamics at 
large and small scales. Further they find that the as- 
pect ratio of local avalanches is consistent with recent ex- 
perimental advances of multiscale roughness analysis Q ■ 
Experimentally, Grob et al. [3l| have, through the termi- 
nology of seismic catalogs, been able to compare the dy- 
namics of interfacial crack propagation to what is found 
in shear rupture for earthquakes. 
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Most of the previous studies mentioned in the above 
paragraph address only rapid event statistics, for a frac- 
ture propagation that is forced by the imposed bound- 
ary conditions (critical fracture propagation). What we 
present here is more elaborate and general in the sense 
that we consider intermittency in both high and low ve- 
locity regimes of crack propagation using two different 
methods of external loading: 1) constant opening veloc- 
ity of the crack and 2) creep relaxation of a crack main- 
tained at a constant opening distance. While it is easy 
to imagine that these different boundary conditions will 
give a very different global behaviour, we are surprised 
to find that the local dynamics is similar in every re- 
spect. This is shown by statistical analysis of high and 
low velocity events, referred to as depinning and pinning 
clusters respectively, and by considering the autocorrela- 
tion of the velocity field. The vanishingly small timecor- 
relations have been related to the time evolution of the 
width of the fracture front [l7|. We see that it follows 
simple diffusion growth. Another important finding is 
that the pinning and depinning size distributions are de- 
scribed by the same power law exponent. Moreover we 
propose a relationship between the different power law 
exponents describing the fracture process, thus linking 
velocity fluctuations with spatial avalanches. 

This paper is organized as follows: In Sec. |TT] we de- 
scribe in detail the experimental setup, including sample 
preparation, loading conditions and optical setup. We 
then present the results in Sec. Mil starting with the 
distribution of local velocities along the fracture front 
(Sec. IIII A|) . In Sec. IIIIBI we obtain the autocorrelation 
functions in time and space for these velocities. Finally in 
Sec. IIII Cl we give the main statistical analysis of spatial 
clusters that we eventually show to be linked to the local 
velocity distribution in Sec. IIII A1 Section HVl summarizes 
the paper with concluding remarks. 



II. EXPERIMENTAL SETUP 
A. Sample preparation 

The experimental setup d, 0, [3l| is shown in Fig. [Tj 
The fracture sample is made out of two transparent Plex- 
iglas (PMMA) plates: a thicker plate with dimensions 
(Zi, w\,hi) = (30, 14, 1) cm and a thinner plate with di- 
mensions {h,W2,h 2 ) = (30,10,0.4) cm for the length, 
width, and thickness respectively. The plates are then 
sandblasted on one side using glassbeads ranging between 
50 /im and 300 /im in diameter. Sandblasting introduces 
random roughness on the originally "flat" surface. This 
causes light to be scattered in all directions from these 
microstructures, hence transparency of the plate is lost 
and it becomes opaque. The plates are then clamped to- 
gether in a pressure frame, with the sandblasted sides fac- 
ing eachothcr. The pressure frame is made of two parallel 
aluminum plates, exerting a normal homogeneous pres- 
sure on both sides of the PMMA. Finally, the pressure 



frame is put in a ceramic temperature controlled oven 
at 205 °C for 30 — 50min. This annealing or sintering 
procedure creates new polymer chains between the two 
plates and the resulting PMMA block is now fully trans- 
parent. The new layer created between the two plates are 
weaker than the bulk PMMA, so that we obtain a weak 
plane with quenched disorder in which the fracture can 
propagate. This system is ideal for direct visual observa- 
tion since the fractured part of the sample immediately 
becomes opaque whereas the unfractured part remains 
transparent. The sharp and high contrast boundary be- 
tween transparent and opaque parts thus defines the frac- 
ture front. 
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FIG. 1. Sketch of the experimental setup. Two PMMA plates 
with dimensions (Zi, Wi,hi) = (30, 14, 1) cm and (Zi, W2, Z12) = 
(30, 10, 0.4) cm are sintered together, creating a weak in-plane 
layer for the fracture to propagate. Fracture is initiated by 
lowering a cylindrical press bar, controlled by a step motor, 
onto the lower plate. The uncracked part of the sample is 
transparent, whereas the cracked part has lost transparency 
hence creating a good contrast at the fracture front. The 
fracture front is imaged from above by a digital camera. The 
deflection d (z-direction) between the plates is indicated in the 
lower panel. The fracture plane is (x, y), where the a;-direction 
is transverse to the average direction of fracture propagation 
whereas the y-direction is parallel to the average direction of 
fracture propagation. 

The rough surface generated by the sandblasting tech- 
nique depends on the volume flux of the beads, the kinetic 
energy of the beads, the bead size, and the total time of 
the sand blasting. It is important to note that there is no 
obvious direct link between the bead size and the charac- 
teristic size of the disorder. The rough surface will after 
annealing give local toughness fluctuations. The strength 
of these fluctuations will depend on the sintering time. 
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The relationship between the disordered morphology of 
the plates and the toughness fluctuations is very diffi- 
cult to access experimentally. However we know that the 
toughness fluctuations will change when the disorder of 
the plates changes 0- In [H[ a white light interferom- 
etry technique was used to measure the rough surface, 
sandblasted with 50 — 100 (im particles, and it was found 
that the local heterogeneities had a characteristic size 
of ~ 15 fim. Other samples have been studied through 
a microscope @ where the random position of the de- 
faults and the maximum size of the defaults was seen to 
roughly correspond to the bead size ~ 50 /im. However 
we emphasize that the image pixel resolution is smaller 
(~ 1 — 5/im) and the largest length scales considered 
(~ 10 3 /Ltm) is much larger than the sample disorder. 

Two different PMMA samples, characterised by the 
glass bead diameter, have been used in our experiments. 
Sample #1 has been sandblasted with 100 — 200 /im 
beads whereas sample #2 has been sandblasted with 
200 — 300 fim beads. Both samples where sintered in the 
oven for 50 min. 



B. Mechanical setup and loading conditions 

The thick plate of the PMMA block is mounted on a 
rigid aluminum frame, also containing a camera setup 
for imaging. Mode-I fracture is induced by a normal dis- 
placement of the thin plate pushed by a cylindrical press 
bar, as shown in Fig.[T] Indicated is also the definition of 
our coordinate system, where (x, y) is the fracture plane: 
the x-direction is transverse to the average direction of 
fracture propagation whereas the y-dircction is parallel 
to the average direction of fracture propagation. The de- 
flection d is defined as the plate separation at the position 
of the press bar. A bit of glycerol is put on the contact 
between the plate and the press bar to reduce any friction 
and prevent shear loading. The press bar is mounted to 
a force gage on a vertical translation stage controlled by 
a step motor, so that it can be moved up and down in 
the z-direction. Through the force gage we are able to 
monitor the force exerted on the lower plate during an 
experiment. To ensure a homogeneous loading, all com- 
ponents of the experimental setup are mounted on a rigid 
plane aluminum plate and leveled. Particularly, a level 
is used on the thin plate to ensure that it is perfectly 
horizontal. If not, adjustments are made to make it so. 
This is also done with the press bar, thus any gradient 
in the loading should not exist. 

We use two sets of loading conditions: 1) The imposed 
deflection d (see Fig. [1) as a function of time t is given 
by 



d(t) 



v p t 



(1) 



where v p is the velocity of the press bar. Throughout 
the experiment we can measure the force F on the lower 
plate at the position of the press bar. As an example, the 
force development during an experiment is shown in the 



upper panel of Fig. [2ja) . Initially there is a period of lin- 
ear increase, corresponding to pure elastic bending of the 
lower plate. At some point, indicated by the dashed line, 
linear behaviour is deviated and fracturing occurs. Af- 
ter some transient period, the force decays only slowly in 
time as the fracture propagates in the sample. The corre- 
sponding linear increase of the deflection is shown in the 
bottom panel. We will refer to these loading conditions 
as constant velocity boundary conditions (CVBC). 
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FIG. 2. (a) Constant velocity boundary conditions (CVBC). 
Upper panel shows the force development F(t) on the lower 
plate as it is bent by the pressbar. The dashed line indicates 
the onset of fracturing. Lower panel shows the linear increase 
of the deflection d(t). (b) Creep boundary conditions (CBC). 
Same as in (a) but F(t) and d(t) are in semilog scale. The 
short solid line in the upper panel indicates the onset of frac- 
turing, whereas the dashed line indicates the time at which 
the pressbar is stopped and maintained in a constant position 
according to Eq. ([2} . 



2) The deflection is given by 



d{t) 



v p t for t < t 
const, for t > t 



stop 
stop 



(2) 



where t s t op marks the time at which the step motor con- 
trolling the pressbar is switched off, i.e v p = 0. We will 
refer to these loading conditions as creep boundary condi- 
tions (CBC), since it is seen that the fracture front con- 
tinues to propagate at "creepingly slow" velocities over 
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several days after t stop - An example is shown in Fig.[2jb), 
where we see a logarithmic decay of the force while the de- 
flection is maintained constant. Motivated by the differ- 
ent global behaviour of the fracture in CVBC and CBC, 
we have performed experiments using both loading con- 
ditions to study the local dynamics. 



C. Optical setup 



A small central region, at the millimeter scale, of the 
front propagation is followed in time using a high speed 
digital camera mounted on a microscope. The large 
width of the bent PMMA plate (10 cm) ensures that the 
central region of interest is not influenced by finite size 
effects (see Fig. [3]). In one experiment between 12 000 
and 30 000 frames are captured using either the Photron 
Fastcam- Ultima APX (512 x 1024 pixels) or the Pixelink 
Industrial Vision PL-A781 (2200 x 3000 pixels). High- 
resolution images (~ 1 — 5/im/pixel) are captured at 
high frame rate relative to the average propagation ve- 
locity of the crack front (see Table HJ. This is important 
as the local fluctuations in velocity can range over several 
decades. As large amounts of data are accumulated, we 
only have the possibility to follow the fracture front over 
short time windows compared to the long-time global de- 
velopment in the examples shown in Fig. [2] Both in the 
case of CBC and CVBC these time windows are small 
enough so that the average propagation velocity of the 
crack front is considered constant. Also for CBC we did 
several experiments with very different average velocity 
(Fig. [2jb)) during the same loading periods. The span 
of the timewindows will of course vary depending on the 
average velocity, but the y-distance (parallel to direc- 
tion of propagation) covered by the crack front is roughly 
~ 500 fxm in all our experiments. Finally, image capture 
is initiated only after onset of the fracture process. 

The obtained grayscale images of the fracture front 
contain two parts: a dark and a bright region, corre- 
sponding respectively to the uncracked and the cracked 
part of the sample. The gray level distribution of the im- 
age thus presents two distinct peaks. Image analysis is 
performed to obtain the coordinates of the fracture front 
line, h(x,t), separating the two regions. This is done 
by thresholding the grayscale image at the local mini- 
mum of the gray level histogram, between the bright and 
dark peak. We then obtain a black and white image from 
which the front can easily be extracted. We always obtain 
a very good contrast between the cracked and uncracked 
part of the sample; the extracted fronts are very robust 
with respect to perturbations in the chosen threshold. 
For a more detailed description of the front extraction 
and image treatment see jy, |31j . 

Fig. [3] shows an extracted front line h(x,t) superim- 
posed on the corresponding raw image. Its roughness is 
due to local pinning asperities of high toughness, created 
as a result of the sandblasting and annealing procedure 
as explained earlier. Occasionally, on small scales close 




FIG. 3. Fracture frontline h(x,t') at some time t', superim- 
posed on the corresponding raw image. Direction of propaga- 
tion is from top to bottom. System size L in the re-direction 
is indicated. The framed raw image corresponds to a tiny 
central part of the full sample as seen in the upper panel. 



to the pixel resolution, the front shows local overhangs 
and is not always a single valued function of x. How- 
ever the number of overhangs per front and the scale at 
which they occur are small; hence we construct the single 
valued front h(x, t) by keeping only the most advanced 
y— coordinate at the front line for a given x— coordinate. 
Arbitrarily we could also have chosen the least advanced 
y— coordinate. Single valued fronts are constructed in or- 
der to simplify the statistical analysis, which has shown 
not to influence the results. 



III. RESULTS 

The rough fracture front exhibits self-affine scal- 
ing properties (l-flTI. [l3l lh"H together with a complex 
avalanche like motion with very large velocity fluctua- 
tions. Due to the large temporal and spatial variations 
in front velocity it is not straight forward to analyze 
the local dynamics by a simple front subtraction proce- 
dure. Therefore we characterize this complex behaviour 
by measuring the local waiting time fluctuations of the 
crack front during its prop agation, following the proce- 
dure introduced first in fl8|. We compute a so called wait- 
ing time matrix (WTM) [30I432T ] , which is a pinning time 
map with elements w, giving the amount of time the front 
is pinned down or fixed at a particular position (x, y) in 
time step units. As explained in Appendix^ the local 
velocity v at a given position is given as v = a/(wSt). 
Using h(x, t) and the WTM, it is then straight forward 
to obtain the local velocities along a fracture front v{x, t). 
Furthermore, by computing v(x,t) for all time steps, we 
build the spatio-temporal velocity map V t (x,t). The av- 
erage velocity (v) is defined as the average over all ele- 
ments of Vt(x,t), i.e the total average over all fronts. 

Presented below are the results of eight experiments 
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(both CBC and CVBC), spanning a broad average prop- 
agation velocity range, where we have characterized the 
local dynamics. The total duration of an experiment is 
within the range of 4 seconds to 7 hours, whereas the 
average distance of front propagation, is ~ 500 /iin in all 
cases. The details of each experiment can be found in 
Table Q] Additionally we will also compare the present 
data to previous experiments from [lj|- 



A. Distribution of local velocities 

A gray scale map of the waiting time matrix is shown 
for a CBC experiment in Fig. [4] Dark regions correspond 
to a high waiting time and thus a low velocity, and vice 
versa for bright regions. The dark low velocity regions 
are seen to occur as irregularly shaped "lines", separated 
by brighter compact regions referred to as high velocity 
avalanches. The wide span of waiting times shown by 
the colorbar, together with their irregular distribution in 
space, is direct visual confirmation of the complex dy- 
namics found in this system. Furthermore, the visual 
impression of the WTM for a CBC experiment compared 
to a CVBC experiment is identical. The similarity of the 
local dynamics in CBC and CVBC experiments is also 
confirmed in our analysis, as we will return to. 

From the local velocities along all front lines Vt(x,t) 
we can compute the normalized probability density func- 
tion (PDF) P(v). By rescaling every local velocity with 
the average propagation velocity v/(v), we obtain a data 
collapse for all experiments as shown in Fig. [5j In this 
figure the results from all experiments in Table U are put 
on top of previous experiments from [l8| . It was found 
that 



P(v/(v)) cx (v/(v)Y 



for v/(v) > 1 



(3) 



with the exponent rj = 2.55±0.15. It is important to note 
that the PDF P(v), computed here directly from Vt (x, t), 
is exactly the same quantity as the PDF of the local front 
velocity v found by estimating the occurrence number of 
each measured velocity on all the pixels in all the fracture 
front line images, as defined in [18|. The result in Eq. ©, 
primarily obtained for CVBC, is now extended to the 
case of creep experiments. It is indeed very stable over 
the different experiments, considering the wide range of 
average velocities. We emphasize that Fig. [5] provides 
quantitative confirmation on the similarity between the 
local dynamics for CBC and CVBC experiments. 

At this point we divide the velocity distribution in two 
and define: a pinning regime for v/(v) < 1 and a depin- 
ning regime for v/(v) > 1, as indicated in Fig. [5] The 
Fig. [5] inset shows the corresponding PDF of waiting 
times P(w/{w)). Through Eq. (|A2[) the two distributions 
are related by P{v)dv = P(w)dw (cf. Eq. ([35])), giving 
P(w/(w)) cx (w/(w)Y>~ 2 for w/(w) < 1. Note that the 
waiting time distribution decays very fast in the pinning 
regime compared to the depinning regime. 



B. Space and time correlations 

The power law distribution of the local velocities con- 
firms the visual impression of a non trivial local dynamics 
of the fracture process. As mentioned earlier, the front 
propagates through high velocity bursts of different sizes. 
An important question is thus how the local velocities 
along and between different front lines are correlated in 
space and time. 

We define the normalized autocorrelation function 
G(Ax) and G(At) for the local velocities on all front- 
lines v(t, x) in space and time as 



G{Ax) 



G(At) 



((v(x + Ax, t) - (v) x )(v(x, t) - (v) x )) x 



((v(x,t + At) - (v) t )(v{x,t) - (v) t ))t 



(4) 



(5) 



where {v) x and a x is the spatial average and standard de- 
viation respectively at a given time in Vt(x,t), whereas 
(v) t and <T t is the temporal average and standard de- 
viation respectively for a given position in Vt(x,t). The 
outer brackets in Eqs. @ and ([5]) denotes an average over 
all different realizations in time and space respectively, i.e 
over all columns and rows in the Vt matrix. 

In Fig.|6] the spatial correlation function G(Ax) is 
shown for all experiments listed in TablelU It is more 
or less evident that correlation functions obtained from 
the same sample are grouped together, independently of 
the average propagation velocity and loading condition. 
By fitting the data with power law functions with an ex- 
ponential cutoff we get 



G(Ax) oc Ax T * cxp(~Ax/x* 



(0) 



where t x = 0.53 ± 0.12 is the average exponent and x* = 
{92, 131} /jm is the average cutoff or correlation length of 
the local velocities in the x-direction, for sample #1 and 
#2 respectively. The quality of the fits is not perfect, as 
can be seen in Fig. [SJ but they represent each group of 
correlation functions fairly well. It is to be noted that 
extracting well defined correlation lengths is not trivial 
in our data. Other estimators of Eq. (0| arc possible to 
use, e.g the power spectrum method. 

In Fig.|7fa) the time correlation function G(At) is 
shown for all experiments listed in TablelU F° r each 
experiment, functional fits analog to Eq. ^ have been 
made. Using the average value of the power law expo- 
nent Tt « 0.43 and different cutoff correlation times t* , a 
good collapse is obtained. We note also that t* is small; 
typically more than two orders of magnitude smaller than 
the duration of an experiment. The inset shows the scal- 
ing of the correlation time with the average propagation 
velocity 



t* = y*/(v) , 



(7) 
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TABLE I. Parameters of the different experiments, sorted after the average propagation velocity of the front (v): System size 
L (x— direction), image timestep St gives the time delay between the capture of two subsequent images, resolution a gives the 
pixel resolution of an image, displacement type denotes the set of boundary conditions used, and the last column indicates the 
sample number. Sample #1 has been sandblasted with 100 — 200 ^m beads whereas sample #2 has been sandblasted with 
200 - 300 ^m beads. 



{v) (/xm/s), L (fim), St (s), a (^im/pixel), displacement type, sample 
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FIG. 4. Waiting time matrix of a CBC experiment, (v) — 1.36 (im/s. The map results from the extraction of 24 576 front lines 
at a rate of 50 fps. Dark regions correspond to a high waiting time and thus a low velocity, and vice versa for bright regions, 
as shown in the colorbar indicating the amount of time (in seconds) the front has been fixed at a given position. Black pinning 
lines are visible, with bright depinning regions in between. The system size L is indicated. 



where y* w 7 /im. The proportionality constant y* has 
the dimension of a length since the scaling exponent 
equals minus unity. This length scale is on the order of 
the pixel resolution a and also within the disorder limit. 
Hence y* is very small and might be influenced both by 
resolution and disorder effects. For comparison we cal- 
culate G(Ay) directly, i.e. the velocity autocorrelation 
in space along the direction of propagation, defined sim- 
ilar to Eq. @ and shown in Fig. [7|b) . We find no power 
law decay in this case but the drop to zero correlation 
occurs between 10 — 20 /im consistently with y* . Correla- 
tion functions from the same sample are shown in similar 
colors (red - sample #2, blue - sample #1). Within the 
interval {a, 20} /im, where a is the image resolution, the 
sample grouping is not so clear as in the case for G(Ax) 
as shown in the inset, but the same initial trend is ob- 
served. This can be attributed to resolution effects and 
the very small correlation lengths. Thus at the time and 
length scales we are looking at, the local velocities are 
considered uncorrelated in the y— direction. 

Since the local fluctuations control the global advance- 
ment of the crack, it is of interest to consider the evo- 



lution of the width of the fracture front in time. This 
growth process is known to depend on the system corre- 
lations. It has been shown previously lfj that uncorre- 
lated growth processes such as simple diffusion, Brown- 
ian motion, etc, can be described by a growth exponent 
a = 1/2. For the present case we define the root-mean- 
square (RMS) value of the front width Ah(t) as 



(Ah(tf)i 
[(h(x,t 



t Q ) - h) - (h(x,t ) - h )]' 



x,t 



(8) 



where h(x, to) is an initial front line and h indicates a po- 
sitional average height at a given time. This differs some- 
what from the usual situation of a front growth from an 
initially flat front. In our case the front width is defined 
as the fluctuations from an initially rough line which cor- 
responds to the geometry of the front at the onset of the 
experiment. The front width is related to the autocor- 
relation of local velocities in time. By rewriting Eq. ([5]) 
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FIG. 5. Distribution of local velocities P(v/(v)) rescaled by 
the average propagation velocity for various experimental con- 
ditions: A range of roughly four decades in average crack front 
velocity including both CBC and CVBC experiments. Sym- 
bols colored red are results from pl|]. A fit to all the data 
for v > (v) shows power law behaviour with an exponent 
—2.55. Inset shows the corresponding waiting time distribu- 
tion P(w/{w)) The exponent transforms in this case to 0.55. 
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FIG. 6. Space correlation functions G(Ax). Functions from 
the same sample are grouped together (sample #2 - filled 
markers, sample #1 - open markers). A power law with ex- 
ponential cutoff has been fitted to each group of correlation 
functions, as indicated by the solid and dashed line for sample 
#1 and #2 respectively (see text). 
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FIG. 7. (a) Time correlation functions collapsed onto ea- 
chother according to a power law with an exponential cutoff 
G(At) = AAt~ Tt exp(-Af/f*). The exponent is n w 0.43. 
Inset shows the scaling between the crossover correlation time 
and average propagation velocity t* ~ 7 /J,m/(v). (b) Space 
correlation function G(Ay) with logarithmic Ay— axis. Con- 
sistently with (a) and Eq. (O, the local velocities become 
uncorrelated after only a short distance (~ 10 — 20 (im) in 
the y— direction. Correlation functions from experiments per- 
formed on sample #2 and #1 have filled and open markers 
respectively. To some extent we see also here grouping of ex- 
periments from the same sample. The difference is however 
not as clear as for the spatial correlations along the trans- 
verse ie— axis (subparalell to the fronts), on the inset showing 
G(Ax) with logarithmic a;— axis. The reason might be that 
the drop to zero correlation occurs close to the resolution scale 
for G(Ay). 



and using that h(t + to) — h(to) = J t * + *° v(t')dt' we obtai 



(Ah(t) 2 } = (\h(x,t + to) - h(x,t )} 2 ) - (t(v)Y 



rt+to rt+t 



to •> to 



(v(n) ■ v(m))dmdn — (t(v)) . (9) 



By substituting n + At = m and using Eq. ([5]) wc get 

rt+to rt+to— n 



(Ah(tY) 



(v(n) -v{n + At))dAt dn .. 



- (*<«» 



to Jto—n 
2 



/ + /<. 



dn 



t+to — n 



dAtG(At) 



(10) 



t —n 



8 



As argued above, we consider the local velocities uncor- 
relatcd in time. The regime where G(At) behaves as a 
power law is very short, and should only affect Ah(t) 
on very small time scales. Thus we approximate the 
autocorrelation function with the Dirac delta function 
G(At) w 6 (At) which gives 



(Ah(t) 2 ) oc t 



(Ah(t) 2 )? ~r 



(11) 



with the growth exponent a = 1/2. Figure [5] shows 
the scaling of the front width as a function of time 
for all experiments. We find indeed a growth exponent 
a = 0.55 ± 0.08 consistent with Eq. (fTTj) . as indicated by 
the fitted dashed line. The large scale crossover is an ef- 
fect of a limited system size in the direction of crack prop- 
agation. Our direct measurement of the growth exponent 
also agrees with the indirect measures in [13, [35| , where 
the front width power spectrum was analysed at different 
times and interpreted in terms of a Family- Vicsek scal- 
ing, with a dynamic exponent k = 1.2 and a roughness 
exponent S = 0.6 giving a = 8 / k = 0.5. 
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FIG. 8. Scaling of the front width as a function of time, 
rescaled with the average velocity. The dashed line corre- 
sponds to (Ah(t) 2 )? oct - 55 . 

Due to the one-to-one correspondence between velocity 
and waiting time [Eq. (|A2[> ] . the above analysis of corre- 
lations could just as well have been performed using the 
latter quantity. Calculating G(Ax), G(At) and G(Ay) 
using w, we obtain approximately the same trends and 
correlation lengths as for v. We turn now to the statistics 
of the dynamical avalanches in the pinning and depinning 
regimes. 



C. Cluster analysis 

1. Spatial map of clusters 

As discussed earlier the local dynamics of the fracture 
front is a mix of pinning lines where the front is fixed or 



only moves slowly, and sudden propagation in high ve- 
locity jumps or bursts. The statistics in both the pinning 
and depinning regimes will be shown to be scale invari- 
ant and characterized by equal scaling exponents. In 
order to study both these regimes we apply a threshold- 
ing procedure to the velocity matrix V(x, y) and obtain 
a thrcsholded binary matrix Vc'- 



V c 



1 for v > C (v) 
for v < C (v) 



for the depinning regime and 



1 for v < ^ (v) 
for v > (v) 



(12) 



(13) 



for the pinning regime. Here C is a threshold constant 
of the orders of a few unities. An example of a thrcsh- 
olded matrix Vc in both regimes is shown, in Fig. |H1 
The geometrical characteristics of the two regimes can 
be seen quite clearly. Depinning clusters (high velocity 
regions) are compact and extend somewhat longer in the 
x— direction than in the y— direction. Pinning clusters 
(low velocity regions) have also a long x— direction ex- 
tension, but are very narrow in the y— direction on the 
other hand. Thus they can be described almost like irreg- 
ularly curved lines in the fracture plane. From Eq. (|12[) 
it is clear that the cluster size decreases with increas- 
ing values of the threshold parameter C in both regimes. 
Obviously one must choose reasonable values of C in the 
two regimes as the number of clusters goes to one and 
zero when C is very small or very large respectively. 

In order for the thresholding of the velocity matrix to 
be consistent, it is important to note that the average 
velocity must be constant in time to avoid clusters from 
being affected by a size gradient. Thus wc ensure that the 
duration of image capture is short enough for the global 
development of the average velocity to be approximated 
as constant for CBC and CVBC experiments. 



2. Size distribution of clusters 

We will denote the size/area of a cluster, for both 
pinning and depinning, S. Figure [TU] shows for C = 3 
the normalized probability density function (PDF) of the 
sizes P(S) respectively for all experiments. There are 
several aspects to emphasize about these figures. First of 
all, the distributions show a power law decay, with a cut- 
off for large sizes S. Furthermore the distributions fall 
on top of eachothcr, meaning that they span the same 
range of cluster sizes, independently of the average prop- 
agation velocity. There is neither no clear indication that 
the PDF cutoffs depend on the correlation length x* . It 
is thus reasonable to average cluster data from all the 
experiments to improve in particular the tail of the dis- 
tribution. Finally, the distributions from both CBC and 
CVBC experiments cannot be distinguished. Thus the 
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FIG. 9. Thresholded matrix Vc (5600 x 1400)/im in the depinning (left) and pinning (right) regime for a CBC experiment with 
(v) — 1.36 jim/a. White clusters correspond to velocities C times larger than (v) for the depinning case, whereas white clusters 
or lines correspond to velocities C times less than (w) for the pinning case. 



distributions seem to indicate that the local dynamics 
are very similar in the two cases, despite very different 
boundary conditions. We will in the following quantify 
the properties of these distributions. 
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FIG. 10. Probability distribution function P(S) for all ex- 
periments using a threshold C = 3. A distribution averaged 
over all experimental conditions is also included for the depin- 
ning (dashed line and circular markers) and pinning regime 
(dashed line and square markers). The pinning size distribu- 
tions have been shifted along the y— axis to enhance visual 
clarity. 

Figure [Til shows the averaged P(S) distributions for a 
threshold range C = 2 — 12 in the pinning regime. It 



is clear that the distributions follow a power law with 
an exponential like cutoff. Furthermore it is evident and 
to be expected that the size of the largest clusters, i.e. 
the cutoff cluster size, decreases with increasing values 
of the threshold level. A similar behaviour is found for 
the PDFs of cluster sizes in the depinning regime, but 
the cutoff size is generally larger due to the cluster ge- 
ometry. In contrast to what was done in [l8||, where 
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FIG. 11. Distributions of pinning clusters P(S) averaged over 
all different experimental conditions, for a threshold range 
C = 2 — 12. Solid lines show the fits corresponding to a 
power law with an exponential cutoff. 
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the distributions were rescalcd by the average cluster size 
(P(S/{S))), we choose to fit the distributions according 
to the function 

P(S) oc S-~< exp(-S/S*) , (14) 

where S* is the cutoff cluster size and 7 the power law 
exponent. This is shown for the pinning regime in Fig.HTl 
where fitted solid lines are plotted on top of the averaged 
experimental data (similar fits have been obtained for 
the depinning regime). We find that in both regimes, 
the cluster size PDF scales with an average exponent 
7 = 1.56 ± 0.04. Using this exponent, and the fitted 
values for the cutoff cluster size we obtain a data collapse 
in both velocity regimes for the full range of available 
threshold values, as shown in Fig. Q21 Furthermore we 




S/S* 

FIG. 12. Collapsed P(S) distributions averaged over all dif- 
ferent experimental conditions for both depinning (upper set 
of data) and pinning (lower set of data). The pinning dis- 
tributions have been shifted for visual clarity. Depinning and 
pinning thresholds are in the range C — 2 — 22 and C — 2—12 
respectively. The dashed and the solid line both have the 
slope 7 = 1.56. Inset shows the scaling between the cutoff S* 
and the threshold C for the depinning (solid line ad = 1.77) 
and pinning regime (dashed line a p = 2.81). 

find a scaling relation between the cutoff cluster size S* 
and the threshold level C, as shown in the inset of Fig. [121 
For the depinning regime it is given by 

S* oc Cr"* , (15) 

where ad = 1-77 ± 0.16. Similarly, we obtain for the 
pinning regime. 

S* oc C~ a » , (16) 

where a p = 2.81 ± 0.23. 

The exponent 7 = 1.56 is somewhat lower but 
consistent with the previously reported value in [l8T | 
(7 = 1.7 ± 0.1), in which the distributions were rescalcd 
by the average cluster size in lack of a pronounced cutoff 
size. A later check using the rescaling as explained in the 



above paragraph does show to lower the exponent also for 
the old data. We would like to mention that our exper- 
imentally obtained exponents 7 and ad are in excellent 
agreement with the recent numerical study of high veloc- 
ity clusters in planar crack front propagation by Laurson 
et al. l 33|. They use an empirical value of ad = 1.8 to 
describe the relationship between the cutoff size and the 
threshold. Their value of the size exponent 7 = 1.5 is ex- 
plained theoretically from the decomposition of a global 
avalanche (collective movement of the front as a whole) 
into local clusters. The experimental equivalent to the 
suggested numerical approach is to study how the fluc- 
tuations of the spatially averaged instantaneous velocity 
(f|(a;, t)) x relates to the distribution of local clusters that 
we observe here. We do not consider global avalanches 
in this study but it is certainly available in our data and 
is a work in progress. 

3. Scaling relations 

The collapse in Fig. [12] shows that the scaling in 
Eqs. (fT5j) and (fl~6| are well satisfied. If we first consider 
the depinning regime, it is possible to relate the expo- 
nents ad and 7 of the cluster size distribution [Eq. (|14j) ] 
to the exponent r\ characterising the spatio-temporal dis- 
tribution of local velocities [Eq. ([3])] . The latter distri- 
bution is obtained from Vt(x, i), i.e the velocity map in 
space and time of all front lines [Eq. (|A3|) and Fig. 1181b)]. 
thus the space-time fraction covered by local velocities 
between v and v + dv is P(v)dv. One may also define 
the spatial distribution of local velocities, obtained from 
the spatial map of local velocities V(x, y) [Eqs. (|A1|) and 
(|A2[)]. denoted R(v). The fraction of (x, y) space covered 
by local velocities between v and v + dv is then R(v)dv. 
As shown in Appcndix|B] there is a relationship between 
these two probability density functions. Using Eq. (|B5[) 
gives 

R( v ) = v ~ v~ ri+1 , for O (v) . (17) 

(v) 

The cumulative distribution of R(v), from a given thresh- 
old C and up to the highest velocity, equals the area 
fraction that these velocities occupy out of the total area 
swept by the fracture front. In terms of threshold level 
we then get 

/>oo 

R c (v >C)= R(v)dv ~ C-"+ 2 . (18) 
Jc 

The same area fraction can also be expressed through the 
cluster size distribution, hence we obtain 

R c (v > C) = ^ (19) 

p 00 

oc N SP{S)dS , (20) 
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where A XtV is the total area in the (x, y) plane where 
the fracture has propagated, TV is the total number of 
clusters, (S) is the average cluster size, and Si ow is 
the pixel size or some other lower cutoff. Substituting 
P(S) = B 5~ 7 exp(— SyS*), where B is the normaliza- 
tion factor, in the above integral, we obtain for (S), 



1 f°° 

-= S-~< exp(-S/S*)dS 
B Js lom 



(21) 



(S)=B S 1 - 7 exp(-5/5*)d5 , (22) 

where S* is the cutoff cluster size. Considering the nor- 
malization factor, we get by substituting x = S/S* 



1 0*1—7 

B ~ b 



x 7 exp(— x)dx . (23) 



Since the lower limit is very small and 7 = 1.56 > 1, the 
power law part of the integrand will dominate and the 
contribution from the upper cutoff is negligible. Thus we 
approximate 



1 



- ra S* 1 ^ / x^dx ~ S* 1 - 1 ^- 

B Js loul /s> 



1-7 

g*l — 7 low ' 



(24) 

which is independent of S*. For the average cluster size 
we then obtain 



(S) oc S* 2 - 1 



A-i 



exp(— x)dx . (25) 



Since 7— 1 = 0.56 < 1, this integral will converge at 
the lower end, to a value independent of Si ow as long as 
Siow/S* < 1. Thus from Eq. (fj"5j). we obtain: 



{S} ocS* 2 ~ 7 (x C-^ 2 -^ 



(26) 



where <7<j(2 — 7) = 0.79. Equation ([2"rJ[) is experimentally 
verified for C > 3, as shown in Fig. [13] 

The number of clusters TV depends on the threshold 
level in a non-trivial manner. This is shown in the inset of 
Fig. H31 We see however that in the interval 3 < C < 16 
the number of clusters can be approximated by 



TV(C) ~ C x 



(27) 



where \ = 0.28. Inserting Eqs. ([T8 jl .([26 ]) and ([27} into 
Eq. (|19[) we obtain the following scaling relation 



^ Q—<Td(2—7)+X 



(28) 



leading to a quantitative link between the exponent of 
local velocity distribution and the exponent of the event 
size distribution: 



77 = (T d (2-7)-x + 2 



Inserting numbers in the above equation (77 = 2.55, 7 = 
1.56, x = 0.28) we get that <Td = 1.88, in good agreement 
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FIG. 13. Average cluster size (S) obtained from the image 
analysis vs. threshold level C. The dashed line shows a power 
law fit for 3 < C < 30, with the exponent — 0^(2 — 7) = —0.75. 
Inset shows the number of clusters as a function of threshold 
level for the various experiments. The dashed lines all have 
the average slope x — 0.28. 



with the empirically found value of ad = 1.8. Strictly 
speaking this result is only valid for 3 < C < 16. If 
we now turn to the pinning regime, we note that from 
Eqs. (fl"5[) and (fl6|) . a p fts aa + 1, allthough we can not de- 
rive it from a theoretical argument. The pinning thresh- 
old values spans a velocity interval {v/(v} < 0.5), in 
which the P(v/(v)) distribution does not follow a power 
law (Fig. [5]) . Thus a similar scaling argument to the de- 
pinning regime, based on simple power law behaviours of 
all dependent variables, is not very likely to hold. 



4- Cluster morphology 

A depinning cluster of size S can be further decom- 
posed into two extension lengths l x - transverse to the 
average direction of front propagation and l y - parallel 
to the average direction of front propagation, by fitting 
a bounding box. A bounding box is the smallest rectan- 
gle that can enclose the cluster, with sides l x and l y as 
shown in the left panel of Fig. [TJ](a). As mentioned ear- 
lier, the pinning cluster geometry can be characterized as 
an irregularly curved line with a much larger extension 
in the x— direction compared to the y— direction. Due to 
this feature, l y is not a good measure, and badly over- 
estimates the y— direction extension. This is shown in 
Fig|T3l(b) where bounding boxes for both pinning and de- 
pinning clusters are shown. Thus for pinning clusters we 
use l x in the x— direction and the average cross sectional 
width l yw as a measure of the y— direction extension, as 
shown in the right panel of Fig. [T?](a). Analysis shows 
that for a cluster of size S, either depinning or pinning, 
the extension lengths have well defined means l Xl l yj and 
(29) lyw increasing monotonically with 5. Note here that the 



bar denote the mean only over a narrow range of S and 
is not the overall mean. The corresponding standard de- 
viations are small and proportional to these means. Due 
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ered equal in both velocity regimes, a y — 0.41 ± 0.06 
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FIG. 14. (a) Left panel shows a bounding box with sides 
l x and i H embedding a depinning cluster. In this case the 
bounding box is a good measure of the linear extension of 
the cluster. Right panel shows a bounding box embedding a 
pinning cluster. In this case l x gives a reasonable linear extent 
measure, however l y does not, due to the irregular curvature 
(somewhat exaggerated in the figure) and to the narrow width 
in the y— direction. To characterize this width we use instead 
the average cross sectional width l yw . (b) Upper and lower 
panel show bounding boxes for depinning and pinning clusters 
respectively from one experiment. 



to the different definitions of l y and l yw , their absolute 
value cannot be compared directly. From analysis we 
find that, after an initial transient, l y and l yw do scale 
similarly but with different prefactors for depinning clus- 
ters. This is a consistency check between using either a 
bounding box or the cross sectional width to describe the 
y— direction extension. Thus l yw is a reasonable measure 
for the y— direction extension of pinning clusters. 

Figure [15] shows the scaling of the different extension 
lengths with the cluster size in the two regimes. In all 
cases there are differences between small (pixel resolu- 
tion up to S ~ 100/im 2 ) and large scale behaviour. In the 
case of depinning, for small S values, l x and l y scale more 
or less similarly indicating that clusters are isotropic at 
these scales. In the case of pinning, l yw is very small and 
stays constant while l x scales almost like the depinning 
cluster size. This is consistent with the characteristic lin- 
ear geometry observed in the pinning regime. However 
the small scale behaviour ranges only over one decade, 
and might be affected both by resolution and disorder ef- 
fects, so we do not have much information at these scales. 
The large scale behaviour spans close to three decades in 
S and displays robust scaling in all cases. From Fig. [15] 
we obtain the following relationship between extension 
lengths and cluster size 

l x oc S a * , l y oc S a « , l yw oc S a *™ (30) 

for S > 100 fim 2 where a x = 0.62 ± 0.04 is consid- 



in the depinning regime, and a y 



0.34 ± 0.05 in the 



depinning regime. The exponents in both regimes con- 
firm the visually observed anisotropy of cluster extension. 
Note also the very small y— direction maximum exten- 
sion (l yw ~ 25 /mi) of pinning clusters, resulting from a 
small proportionality factor in the scaling relation. Fur- 
thermore we obtain approximately from the exponents in 
Eq. ([30]) that S ~ l x l y ~ l x lyw, meaning that the ratio of 
the approximated area from the extension lengths to the 
real cluster area is scale independent. From Eq. ([30]) we 
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FIG. 15. Linear extent of pinning and depinning clusters as 
a function of cluster size for the full span of threshold levels 
and averaged over all experimental conditions. The slopes 
of the different fitted lines (dashed - pinning clusters, solid - 
depinning clusters) are indicated in the caption. Note that 
there in all cases are initial transients up to S ~ 100 fim 2 . 

get the following x— and y— direction aspect ratio: 

K ^x v ^ x ' lyw l" vw x , (31) 

where a y /a x = 0.66 and a yw /a x = 0.55 for the depin- 
ningand pinning regime respectively. It was suggested 
in [181 ] and in |36| that a y /a x could be another mea- 
sure of the roughness of the self-afHne fracture front, in 
agreement with previous experimental measurements of 
the roughness exponent. However, in a very recent ex- 
perimental work Q on planar crack growth, there has 
been two roughness exponents observed acting at differ- 
ent scales; a smallscale roughness with exponent ~ 0.6 
and a largescale roughness with exponent ~ 0.4, with a 
crossover depending on the fracture toughness fluctua- 
tions and the stress intensity factor. This trend has also 
been seen for the aspect ratio of depinning clusters in 
the simulation study by Laurson et al. (33j . In the ex- 
perimental case on the other hand, considering that the 
length scale of this roughness crossover are comparable 
with the l x range in our case, we find no traces of such 
behaviour in the aspect ratio of depinning clusters. This 
point thus warrants further consideration. 

Finally, we discuss the marginal distributions of the 
extension lengths, i.e. for all cluster sizes, in the two 
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FIG. 16. (a) Collapsed P(l x ) distributions averaged over all 
different experimental conditions for both depinning (upper 
set of data) and pinning (lower set of data). The pinning 
distributions have been shifted for visual clarity. Depinning 
and pinning thresholds are in the range C — 2 — 30 and C = 
2—12 respectively. The solid and dashed lines both have 
the slope j3 x — 1.93. Inset shows for the case of depinning 
the threshold dependence for the unsealed distributions, (b) 
Collapsed P(l y ) and P(l yw ) distributions averaged over all 
different experimental conditions for the depinning (upper set 
of data) and pinning (lower set of data) regime respectively. 
The pinning distributions have been shifted for visual clarity. 
Thresholds are in the range C = 2 — 30 and C = 2 — 12 for 
depinning and pinning respectively. The solid line has the 
slope f3 y = 2.36. Inset shows for the case of depinning the 
threshold dependence for the unsealed distributions. 



regimes denoted P(l x ), P(ly), and P(l yw )- For clarity 
we mention again that l x scales similarly with S in the 
two regimes only separated by a small difference in the 
proportionality factor, whereas l y describing the depin- 
ning regime, and l yw describing the pinning regime, are 
treated separately. The insets in Fig. [TBT a) (b) show the 
extension length distributions P(l x ) and P(l y ) respec- 
tively in the pinning regime. The corresponding pin- 
ning cluster distributions display similar behaviour, ex- 



cept that the P{l yw ) distribution is entirely dominated 
by a cutoff function. This is due to the very narrow 
y— direction span of pinning clusters. We define the fol- 
lowing distributions for the extension lengths 

P(l x ) cx l-?*D(l x /l*) (32) 

P(ly) CX l-^D(ly/l* y ) (33) 

P(ly W ) oc ly^D(l y /l y ) , (34) 

where D(x) is some cutoff function decaying faster to 
zero than any power of l x , l y or l yw when x > 1 and con- 
stant otherwise. The (3 exponents above can be predicted 
from our previous results for the cluster size distribu- 
tion. From statistics we know that the relation between 
the PDFs of two random variables b and c, one-to-one 
related, can be expressed as 

P(6) = P(c)J. (35) 

In our case S, l x , l y , and l yw is not one-to-one related, 
but since the means l x , l y and l yw have only small stan- 
dard deviations, the PDFs P(l x ), P(ly), P(l yw ) should 
at least be approximated by Eq. (f3"5j) . For P(l x ) we get 
by inserting Eq. flU) and Eq. into Eq. ([35]) 

Px = , (36) 

a x 

where /3 X = 1.93. Similarly we obtain /3 y = 2.36 and 
f3 yw = 2.65. For the depinning regime we obtain for the 
cutoffs in Eqs. (02) and (J33j by using Eqs. (H) and ([30): 

I* cx C- ad > a * , l*y cx Cr""* . (37) 

For the pinning regime we obtain for the cutoffs in 
Eqs. ((32) and ([34]) by using Eqs. (fTS) and (|5Q )l : 

l* x cx C- a ^ , l* yw cx C- a " a ^ . (38) 

The extension length distributions in both velocity 
regimes are collapsed according to Eqs. (|32M34[) as shown 
in Fig. [T6l a)(b). In the x— direction, transverse to the 
direction of crack propagation, the distribution in both 
regimes scale with the same exponent, similarly to what 
was found for the cluster size distribution. The only dif- 
ference between the two distributions is the proportional- 
ity factor in the cutoff length, as explained earlier. We see 
that along the direction of crack propagation the depin- 
ning [P(l y )] and pinning [P(l yw )] distribution are quite 
different, in the sense that all power law behaviour is sup- 
pressed by the cutoff function in the latter distribution. 
This is understandable since the span of l yw values is no 
more than one decade. 

In Sec. IIII Bl we discussed various correlation functions 
of the spatio-temporal velocity field. In particular it was 
seen that the local velocities had correlation lengths of 
the order ~ 100 (im and 10 /im in the x— and y— direction 
respectively. One would expect the correlation lengths in 
some sense to control the extent of pinning and depinning 
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clusters. This dependence is non trivial since a cluster 
in this context is artificially constructed by threshold- 
ing the velocity field. No clear relation is found between 
the cutoff size of the pinning and depinning clusters, and 
the correlation length extracted from the autocorrelation 
function of the velocity field. However, since the clusters 
are obtained from thresholded velocities, it is also pos- 
sible to look at the autocorrelation function of thresh- 
olded velocities, rather than the one of the continuous 
velocity signal. In ongoing work we consider such cor- 
relation functions Gc(A.t) [Eq. (J2)], obtained from dis- 
cretized signals vc(x,t) where the local velocities along 
each front line are now thresholded with a threshold C 
according to Eq. (|12|). Preliminary analysis indicate the 
existence of a correlation length roughly proportional to 
I* [Eq. (|37[) ]. meaning that both quantities evolve simi- 
larly with the threshold C. 

Furthermore, in the x— direction we could see clear 
sample differences in the correlation lengths, even though 
they were within the same order of magnitude (Fig. [5]). 
Analysing carefully both size and extension length distri- 
butions of individual experiments, and not average distri- 
butions as presented above, we could not recognize such 
trends. In this respect it is also important to mention 
that for individual experiments, the cutoff behaviour in 
the distributions are not well pronounced due to the lack 
of large scale statistics. Even when considering the above 
limitations, we can say that the geometry of pinning lines 
are qualitatively consistent with the observed correlation 
lengths. Thus it seems that the vanishingly small corre- 
lation length in the y— direction, describes the low value 
part of the local velocity distribution. 



IV. CONCLUSION 

The local dynamics of an in-plane mode-I fracture have 
been studied experimentally using high resolution moni- 
toring of the front line advances. Indeed the transparency 
of the PMMA enable us to follow the fracture process us- 
ing a high-speed camera. Fracture is induced by fixing 
the upper plate, while applying a force on the lower plate 
from a press bar controlled by a step motor. Experiments 
are performed using two sets of boundary conditions: 1) 
constant driving velocity on the pressbar, giving a lin- 
ear deflection in time between the plates (CVBC) and 2) 
fixed deflection between the plates (CBC), resulting in a 
slow creep motion of the fracture front. 

Disorder is introduced in the fracture plane by a sand- 
blasting and sintering procedure, resulting in heteroge- 
neous fluctuations of the local toughness. The compe- 
tition between the toughness fluctuations and the long 
range damping elastic forces results in a rough fracture 
front with self affine scaling properties. In this study 
we have considered the local dynamics of the fracture 
front over a wide range of average propagation veloci- 
ties (0.028 < (v) < 141)^m/s. The local velocity field 
is obtained through the waiting time matrix and gives 



a spatio-temporal distribution with a large power law 
tail for high velocities described by an exponent — r\ = 
—2.55. The fracture front advance, displays pinning and 
avalanches with a broad range of velocity scales. Our 
results show that the local dynamics is similar in every 
respect for the two different boundary conditions. This is 
an important and non-trivial result considering the very 
different behaviour in the global large scale propagation. 
Additionally, no dependence on the average propagation 
velocity for different experiments is found. 

The average autocorrelation of local velocities have 
been studied in both spatial directions, and also in time 
along the direction of crack propagation. We find that 
the velocities are correlated up to ~ 100 /xm transverse 
to the direction of crack propagation, and ~ 10 /urn, i.e. 
close to the spatial resolution, and thus uncorrelated in 
the direction of crack propagation. Within these gen- 
eral trends we have seen that there are differences in the 
autocorrelation function from sample to sample, but no 
dependence on the loading condition or average propa- 
gation velocity. Relating the autocorrelation of velocities 
in time to the evolution of the front width gives a growth 
exponent of a = 1/2 similar to simple diffusion, a process 
such as Brownian motion. 

The local dynamics have been studied through a statis- 
tical analysis of local avalanche events. We have observed 
that the cluster properties are independent of both load- 
ing conditions and average velocity of the crack front. 
The depinning cluster size distribution show scale invari- 
ance, described by an exponent —7 = —1.56, in agree- 
ment with previous experimental [18j and numerical re- 
sults [HI, [33 ■ Surprisingly the same result is found also 
for the pinning regime. Furthermore, we have in this 
study seen that the cluster size distribution scaling is 
truncated by an upper cutoff, depending on the threshold 
value. We have shown that the cutoff essentially is con- 
trolled by the total distribution of local velocities. Par- 
ticularly for the depinning regime we have obtained a 
scaling law relating the cluster size exponent 7 to the 
exponent 77 describing the local velocity distribution. 

Clusters have in both velocity regimes been fur- 
ther decomposed into extension lengths in the x— and 
y— direction. We have demonstrated that the distribu- 
tions of these extension lengths are consistent with their 
size distribution. The aspect ratio of depinning clusters 
follows a power law with the exponent a y /a x = 0.66 in- 
dicating that the clusters are anisotropic and extending 
longer transverse to the direction of propagation than in 
the direction of crack propagation. We have yet to ob- 
tain experimentally a relationship between the extension 
of depinning clusters and the roughness of the fracture 
front. This is a topic that warrants further work. 

The pinning clusters were found to display a very 
strong anisotropy, extending far in the x— direction as 
opposed to the very short y— direction extension. This is 
qualitatively in agreement with the found velocity corre- 
lation lengths in the two directions, thus indicating that 
these lengths describe the spatial correlations of low ve- 
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locities. 



Appendix A: The Waiting Time Matrix 

The waiting time matrix (WTM) is a robust procedure 
that enables a comparison of both different experiments 
at different time and space resolution, and also with nu- 
merical simulations of similar systems. It can be applied 
to any propagating interface |3JJ-|32|, and is particularly 
suited for estimating the local velocity of pinned inter- 
faces which are dominated by low speeds. Below, we will 
explain the procedure in detail. 

The coordinates of the extracted front lines h(x, t), in- 
troduced in Fig. [3j can be represented in matrix form 
as: H (x, h{x, t)) = 1 and elsewhere, with a matrix size 
equal to the captured image size. We define the WTM 
W as the sum of all front matrices H, 



W(x,y)=Y,H(x,h(x,t)) 



(Al) 



where the sum runs over all discrete times t. Note that W 
is an integer matrix, so to get the true waiting time, the 
time step St must be multiplied to each matrix element 
w. An example of front line addition is shown in Fig. 1171 



0000000000000 
1000000100001 
0100101010010 
0011010001100 
0000000000000 
0000000000000 

0000000000000 
2000000100001 
0200101120021 
0022111002200 
0000010000000 
0000000000000 



t 3 



0000000000000 
2000000100001 
1200101120032 
0133112112 3 00 
00001 2 0001000 
0000000000000 



FIG. 17. Example of the computation of the waiting time 
matrix W(x,y) [Eq. (|A1|) ]. All fronts are added to an origi- 
nally empty matrix in time step unit. Indicated above is the 
addition of front lines in three timesteps ti (red), ti (blue), 
and tz (green). 

From above it is clear that the WTM procedure gives 
a spatial map that accounts for the amount of time spent 



by the front at a given pixel, thus reflecting the local dy- 
namics of the interface. However, avoiding holes in the 
WTM, implies a high enough sampling rate, so that the 
movement of the front position is at maximum one pixel 
between two subsequent images. Second, it also requires 
a small noise from the imaging device. Finally, care must 
be taken in preparing the sample. Indeed, impurities and 
surface scratches are not transparent but rather reflect 
light and may thus artificially alter the extracted front 
shape. In our case, experiments are devised so that the 
front is propagating in a steady manner both before and 
after the short interval of image capture. To avoid tran- 
sient effects at the beginning and at the end of the image 
recording, we typically clip between 200-500 front lines 
in the start and end of the generated WTM. 

From the WTM we can construct the local velocity 
matrix in space V(x, y). Matrix elements represent the 
normal speed of the fracture front at the time it went 
through a particular position 



1 a 

w St 



(A2) 



From the local velocity matrix V(x, y), we can also obtain 
the local velocity along each front h(x,t) 



v{x,t) = V(x,h(x,t)) . 



(A3) 



By computing v(x,t) for every time step, we build the 
spatio-temporal velocity map Vt(x, t). We then define the 
average propagation velocity of the front (v) as the av- 
erage taken over all elements in the matrix Vt(x,t). The 
development of the front in time for a given x— position 
is shown in Fig.[T8](a), also indicating how the velocity 
is approximated from the WTM. One realization of the 
local velocity fluctuations along a front line is shown in 
Fig.QKb). 



Appendix B: Velocity PDF transformation 

In transforming from the spatio-temporal map Vt(x, t) 
[Eq. (|A3)) and Fig. H5Ib)] to the spatial map V(x,y) 
[Eqs. (SD and (E2J] with the PDFs P(v) and R(v) re- 
spectively, we can express the space travelled through 
at speed v over a time dt as dy = v dt. The area in 
(x,y) space where the front travels at speed u between 
v and v + dv corresponds to the total area of fracture 
propagation, A x ^ y , multiplied by the fraction of the area 
corresponding to this speed: 



v<u(x,y)<v+dv 



dxdy = A Xt yR(v)dv 



(Bl) 



This area is related to the area covered by the fronts trav- 
eling at that speed in the spatio-temporal map, expressed 
using the variable change between y and t: 



dxdy 



v<u(x,y)<iv-\-dv 



dxvdt 



<u(x,t)<v-\-dv 



(B2) 



16 



(a) 




126 



1120 1160 1200 1240 1280 1320 
t (# timesteps) 



Eventually, this last area is directly related to the distri- 
bution P{v), with the same argument as for the spatial 
map: denoting A Xjt the total area of the spatio temporal 
map, we can write 



L 



v<u(x,t)<v-\-dv 



dxdt = A Xjt P{v)dv (B3) 



Inserting Eqs. (|BTJ) and (|B3J) into Eq. flB2j) leads to 



(b) 
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A x<y R(v)dv = A x<t P{v)vdv 



(B4) 



Furthermore it can be shown that A XjV /A Xt t = (v), thus 
eventually 



v P(v) dv = (v)R(v) dv 



(B5) 



FIG. 18. (a) Pixel level zoom in of a frontline h(x',t) at 
a given position x as function of time. Indicated are three 
waiting times Wi, u>2 and 11)3 separated by a one-pixel jump. 
As an example, note that all captured fronts from h(x', to) to 
h(x' , to + Wi) is given the same constant velocity Vi oc 1/wi in 
making the jump from pixel 127 to pixel 128 along the y— axis. 
This approximation means that the front position increases 
linearly during this time interval (indicated in red), (b) Local 
velocity fluctuations v(x,t') along the frontline h(x,t'). 
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